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Abstract 

We calculate the entanglement entropy of the SU{N) Yang-Mills gauge theories on the lattice 
under the strong coupling expansion in powers of /3 = 2N/g‘^, where g is the coupling constant. 
Using the replica method, our Lagrangian formalism maintains gauge invariance on the lattice. 
At 0{fH) and 0{/3^), the entanglement entropy is solely contributed by the central plaquettes 
enclosing the conical singularity of the n-sheeted Riemann surface. The area law emerges naturally 
to the highest order 0{j3^) of our calculation. The leading 0{(3) term is negative, which could in 
principle be canceled by taking into account the “cosmological constant” living in interface of the 
two entangled subregions. This unknown cosmological constant resembles the ambiguity of edge 
modes in the Hamiltonian formalism. We further speculate this unknown cosmological constant 
can show up in the entanglement entropy of scalar and spinor field theories as well. Furthermore, it 
could play the role of a counterterm to absorb the ultraviolet divergence of entanglement entropy 
and make entanglement entropy a finite physical quantity. 
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I. INTRODUCTION 


Entanglement entropy is a measure of the level of entanglement between the degrees of 
freedom in two subregions of a physical system. Besides being a fundamental and mysterious 
property of quantum mechanics and quantum held theory, quantum entanglement is of 
practical use. In some systems, entanglement entropy plays the role of an order parameter 
to characterize quantum phase transitions m, while in others it demonstrates the scaling 
behavior mu near the critical point. In held theory, a widely used method for calculating 
the entanglement entropy is the replica trick j5]. This method calculates the trace of the 
reduced density matrix to the n-th power in the path integral formalism, which amounts to 
computing the free energy of the system on a n-sheeted Riemann surface, or equivalently 
on a “cone” with a conical angle of 2n7T. The entanglement entropy is then obtained as a 
response of the free energy to the change of the conical angle at n = 1. This is similar to 
calculating the black hole entropy by evaluating the response of the free energy of helds with 
respect to the dehcit angle in the Euclidean spacetime [7]. 

In previous studies, while the computation of entanglement entropy for the scalar and the 
spinor helds are considered straight forward, it becomes more subttle for the gauge helds. 
It was found by Rabat [7| that the gauge helds in the black hole entropy method yields an 
extra negative contact term compared to their entanglement entropy. This term is due to 
the interaction of the particles with the horizon, and is believed to be related to the choice 
of the boundary condition while removing the tip of the cone due to the black hole in the 
Euclidean space. It was later interpreted by j8] as arising from incorrect treatment of the 
zero modes. 

As for the gauge helds on the lattice, the entanglement entropy had been studied with the 
Hamiltonian IMl] and the Lagrangian HMD approaches. In the Hamiltonian approach 
for the case of the gauge helds on the lattice |M2], one needs to impose the Gauss law 
or the gauge hxing constraints in order to get rid of the unphysical degrees of freedom, 
and the ambiguity arises in determining which subregion the gauge links on the boundary 
belong to. This leads to the difficulty in decomposing the global gauge invariant states 
into the direct products of those living in each subregions. It was proposed j9] that these 
ambiguities might be compensated by edge modes living in the interface of the subregions 
and determined by the transverse electric helds. This echoes the concept that the contact 
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term in [7] arises from the sources on the horizon. In [9] this negative contact term was 
found to arise from the entanglement of the edge modes. Moreover, in order to overcome 
the puzzle of Hilber space decomposition, m employs various choices of electric, magnetic, 
and trivial centers in the operator algebra in the Hamiltonian formalism, among which the 
trivial center is related to some spacial gauge fixing choice, and hence the entanglement 
entropy thus obtained is not gauge invariant. Besides, the electric and magnetic choices also 
result in different tripartite topological entanglement entropy. In view of these ambiguities, 
we explore a complementary approach—the Lagrangian formalism—to shed light on the 
problem from a different perspective. 

The replica method arises from taking the derivative of the trace of n copies of the 
reduce density matrix by n. Since the density matrix can always be expressed in the path 
integral formalism as long as the quantum field theory is local, the replica method is valid 
for the gauge fields at continuum. We start from the replica method at continuum, and 
then discretize the spacetime into a squared lattice. We divide an infinitely large system 
into two semi-infinite subregions by a flat boundary, and decompose the spacetime into a 
direct product of a 2-dimensional cone with a conical angle of 2n7r (or equivalently an n- 
sheeted Riemann surface) and an ordinary Euclidean space transverse to cone. Then we 
use the Wilson gauge action which has the advantage of gauge invariance on a discrete 
lattice. This action sums over the Wilson loops on the plaquettes, including those on the 
n-sheeted Riemann surface and those on the ordinary Euclidean space. In contrast with the 
previous studies where the conical singularity is placed on the lattice site, we use a different 
discretization setup by locating the conical singularity in the center of the plaquette. As a 
result, the branch lines on the n-sheeted Riemann surface cut across the links, and there 
is no lattice site on the cut. (See Fig. [^) This setup yields two types of plaquettes on 
the squared lattice: the central plaquettes encircling the tip of the cone, formed by 4n 
links, and the regular plaquettes formed by 4 links with no singularity inside. When the 
conical angle, or n, changes, only the central plaquettes (i.e. those plaquettes with 4n links) 
respond to this change. As a result, entanglement entropy necessarily involves those central 
plaquettes. The fact that all of the central plaquettes live across the interface between the 
two subregions naturally give rise to the area law, which states that the leading contribution 
to entanglement entropy scales as the area of the interface. 

The connection to the area law can be further demonstrated order by order diagram- 
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matically under the strong coupling expansion. Interestingly, we find the leading term in 
the strong coupling expansion negative. However, symmetries of the action allow a two- 
dimensional cosmological constant living in the interface |18) which could provide a positive 
contribution at an even lower order. We speculate that the freedom to tune this unknown 
two-dimensional cosmological constant corresponds to the ambiguities encountered in the 
Hamiltonian approach. We further speculate that the two-dimensional cosmological constant 
can show up in the entanglement entropy of scalar and spinor field theories as well. They can 
play the role of a counterterm to absorb the ultraviolet divergence of entanglement entropy 
and make entanglement entropy a finite physical quantity. 

This paper is organized as follows. In Sec. H, we briefly review the notion of entanglement 
entropy and the replica method. The entanglement entropy of Yang-Mills fields on the lattice 
under the strong coupling expansion is calculated in Sec. HI, and the cancellation of the 
negative contribution by including a 2 dimensional cosmological constant is discussed in 
Sec. IV. Sec. V compares our result with the previous ones obtained by the Hamiltonian 
methods, and Sec. VI summarizes our study. 


II. ENTANGLEMENT ENTROPY AND THE REPLICA METHOD 


Suppose our system occupies an infinitely large and fiat d+1 dimensional spacetime and 
is divided into two semi-infinite subregions A and B. They are divided by an infinite and 
fiat d — 1 dimensional space-like boundary. The entanglement entropy (EE) of a quantum 
theory between the two subregions is defined by the von Neumann entropy. With some 
simple algebra, it can be re-expressed as: 


*S'ee = -tr[p^lnp^] 


A 

dn 


IntrK] 

n—¥l 


( 1 ) 


where pa = trs[p] is the reduced density matrix by tracing out the degrees of freedom in 
region B. This expression is called the replica method because it involves n copies of pa- 
An elegant path integral formulation to compute the entanglement entropy using the 
replica method was first introduced in j6] (see also [I9]). In this set up, one recalls that pij = 
(i j) and Tr[p] is the partition function calculated in finite temperature field theory 

with appropriate boundary conditions (periodic and anti-periodic boundary conditions for 
bosons and fermions, respective) imposed for fields at Euclidean time r = 0 and 1/T, where 
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T is the temperature. Then Tr[p^] can be computed by doubling the period (by imposing 
appropriate boundary conditions at r = 0 to 2/T). Similarly, Tr[p^] is computed by doubling 
the period (from 0 to 2/T) for region A while region B still has the single period (from 0 
to 1/T) as shown in the left plot of Fig. [^A), which is equivalent to performing the path 
integral on a 2-sheeted Riemann surface in the right plot. One can generalize this set up to 
Ti/pd] for an arbitrary n. There is no restriction on the space partition between A and B. 
The sizes of A, B, and T can be either finite or infinite. 

In this paper, we will just concentrate on the simplest case with the sizes of space and 
(Euclidean) time to be both infinite (i.e. T = 0) and the interface between A and B to be 
a flat infinite plane. In this limit, the n-sheeted Riemann surface has a conical structure as 
shown in Fig. [^B) with the time and longitudinal spacial direction (the direction that is 
perpendicular to the interface) lying on the cone while the space on the interface transverse 
to the cone. 

As a result, tr[p( 4 ] becomes a partition function on the n-sheeted Riemann surface, or, 
in our case, a cone with 2n7r conical angle, normalized by n-copies of the partition function 
on the ordinary Euclidean space Z”: 


which ensures that as n 


1, tr[p) 4 ]=l. The entanglement entropy is then given by 


( 2 ) 


*S'ee = - ^ (In Z„ - n In Zi) 




” [In Zi+, - (1 + e) In 


(3) 


Note that n is taken as an integer in the integral of Z„. After one obtains the analytic 
expression for tr[p”], then n can be analytically extended to non-integers to carry out the 
differentiation at n = 1. 


III. CALCULATION ON AN N-SHEETED LATTICE MANIFOLD 

We now discrete a d -|- 1 dimensional spacetime into a squared lattice. The spacetime is 
decomposed into a direct product of a 1 -|- 1 dimensional n-sheeted lattice (which we call 
the parallel dimensions) with coordinates (xy, r) and a discrete d—1 dimensional transverse 
space labelled by x±. We discretize the n-sheeted Riemann surface in such a way that the 
end point of the cut on each sheet (or, the conical singularity corresponding to the tip of 
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(A) 




B 


■ft 



(For n = 2) 







Figure 1: The n-sheeted Riemann surface in the replica trick. (A) illustrates the equivalent geomet¬ 
ric structures of the n-sheeted manifold in the case of n = 2 arizing from the replica trick, on which 
the partition function Zn is computed. The unshaded and shaded parts label the subregion A and 
B, respectively. The subregion B is traced out in the reduced density matrix p^. The 2 sheets on 
the right are at the same x_\_ coordinates, attaching to each other at the conical singularity which is 
represented by the black dot and the dotted vertical line. In the right figure, it is also demonstrated 
explicitly after discretization a central plaquette with 4n edges encircling the conical singularity, 
and the ordinary plaquettes with 4 links located on a sheet and across two adjacent sheets (see also 
Fig. §. The numbers label the order of the link variables to form the plaquettes. When the sizes 
of space and (Euclidean) time are both infinite (i.e. T = 0), this geometry is equivalent to a cone 
with 2n7r conical angle depicted in (B). 

the cone) locates inside a plaquette rather than sitting on a lattice site. It will be clear 
later that our result is independent of the location of the conical singularity as long as it 
is encircled by a plaquette. For simplicity we choose that the conical singularity sits in the 
center of a plaquette. 

As a result of the discretization, there are two types of plaquettes. The ordinary pla¬ 
quettes contain no conical singularity inside, and are composed of 4 links. They make up 
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all plaquettes in the transverse dimensions and most plaquettes in the parallel dimensions, 
and their plaquette variables are denoted by where k indicates their 

location bring on the fc-th sheet, and V signihes the ordered product of the link variables 
Ui forming the plaquette. On the other hand, the central plaquettes enclose the conical 
singularity on each sheet in the parallel dimensions and are composed of 4n links. Their 
plaquette variables are denoted by Ub{x±) = 'PYIi^bUi, where x± is the location of the 
parallel planes in the transverse dimensions. See Fig. [^for the cartoon of these two kinds 
of plaquettes. 

Recall that the partition function of the lattice gauge theory on a one-sheet manifold is 
given by 


Z = I Pf/exp|-/?^ 


1 — — Re trf/n 
N 


a^O 


VA exp - / d'^x 


4^2 


(4) 


where fd = 2N /and □ labels the location of plaquettes. The plaquette variable Ub is 
the local Wilson loop composed of the ordered product of four gauge links, Ub ='P ri/en 
where V indicates the ordered product and Ui is the link variable representing the gauge 
helds. The action recovers the Yang-Mills action in the continuum limit by setting the lattice 
spacing a —0. 

To construct a lattice gauge held system in a general d + 1 dimensions whose 1-1-1 
dimensions is an n-sheeted manifold, we rewrite the partition function in Eq. Q in terms of 
the ordinary plaquettes and the central ones Ub{x±) , such that the partition function 
reads 


Z„ — 


a^O 


exp 


'"EE Retr[[/^^^ - 1] + ^ Re tr[f/H(x_L) “ 1] [ (5) 


N 


k=l □ 


m=l 

n 1 fir n ^ 

/ ‘*'^4EX 


X_L 


m=l 


'-{ 0 } 




^ k,i=i 


P[ VA^ 


m) 


m=l 


exp 




where we have assumed the boundary condition 


R^^^(x±,X|| = 0) = F(')(xx,X|| = 0), {k ^ /), (6) 
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Central Plaquette Regular Plaquettes 


Figure 2: One central plaquette (left) enclosing the conical singularity v.s. three regular plaquettes 
on three different sheets (right) in the case of n = 3 sheets. In the left figure, the dotted line 
represents the conical singularity to which all three sheets attach, and the wavy lines are cuts on 
each sheet. The central plaquette has 12 edges and it takes Gvr to go around it, while each regular 
plaquette has 4 links. 


i.e. the field strength on each sheet is assumed be equal. This condition is natural if we 
impose 27r rotation symmetry on the n-sheeted surface. We have introduced an extra 1/n 
factor to the central plaquette terms, because the 4n gauge links of the central plaquettes 
encircle the Tqi flux over an area na^, and give rise to a factor oc to the action, 

inconsistent with the contribution from the transverse plaquettes scaling as na^. Therefore 
a factor of 1 /n is introduced to compensate this effect. 

Expanding Eq.([^ to the second order of /3, we have 


Zn — 




vu^ 


ra) 


/S 


m=l 
2 ^ 


1 + 


N 


f/r 


(k) 


+ 


fc=l □ 


A 

nN 


^Re tr [[/□(x_l)] 




2iV2 


Retrf/f 


r(^) 


Re tr A 


k,l=l 


1 


[Re trf/ 0 (xij_)] [Re trf/H(x 2 _L)] + cross-term -|- 0(/9^ 


(7) 


Xl±,X2_L 


Here A/” is the number of plaquettes on a single sheet for plaquettes not encircling the 
conical singularity, including both the parallel and transverse plaquettes. 7 V_l is the number 
of plaquettes encircling the conical singularity. 

























With RetrU = UtrU + trW), we first evaluate the leading contribution of the Wilson 
loops from f VU trf/trf/h 




1 


2 iV 2 2 


EE 

k,l=l ■ 


/32 1 1 


2 iV2 2 n? 


E 


Xl±,X 2 _L 


Yl 

_m=l 

n 


trf/r 


(fc) 




.k=l 


[trt/0(xi_L)] trl7^(x2±) + 0(/?^) 


= e-z^nAt-^A^x \ I + + 0(/33) . 


4Ar2 4Ar2 n2 

Here we have used the Haar measure 


dUUij = 0=1 dUUjj, 


dUU.jUl = 


N 




such that 


J vutiuShiug^^ = 5ku 
VU trt/0(xi_L) trt/^(x2±) = h 


XiX2 • 


Note that 0 and ([T^ are independent of n. Taking the combination 


In — n In Zi = —/9A/j_ 
and then one obtains the Renyi entropy 


1 

- n 

n 


+ 


/?Wj 


Vn 2 


--n +0(/?3), 


Sn = -PAf± 


1 + n /3‘^N‘± f 1 + n + 


+ 


n 




+ 0{(3^ 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


n 4Ar2 

This means no matter how we place the lattice, as long as the conical singularity inside the 
central plaquettes (i.e. the central plaquettes exist), the result will be the same. But if 
one chooses to put the conical singularity on a lattice site with no plaquette encircling the 


conical singularity, then there is no central plaquettes to offer the ^ and ^ factor in (13), 
but the factors n appears instead and effect of lnZ„ is canceled by that from nhiZi. 

Another non-trivial contribution from the Wilson loops is J T>U (tiU)^. Using the iden¬ 
tities 



^nJn 

1 

1 au ^1232 " 

^NjN 

1 


spq-’S 

rag--s 


(15) 

(16) 
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one finds such contribution are of 0{/3^): 


j VUitiUuf = 1, 

= 1 , 


(17) 

(18) 


= 


vu^ 


m) 


m=l 


/S 


<N 


N\N^ 


^J^Retr 


. k=l □ 


f/r 


(fc) 


n 


^Re tr [t/H(x_L)] 


N 


x± 


2 f AO. 

" CaV + ^ + cross-term 


(19) 


which is subleading for N > 3. For SU(2) theory, we have an additional contribution, 


InZf 2) = + ;^nAr + 


For SU(3) theory, it is 

,.rZf-) = + (21) 

Putting together (13) and (19), the Renyi entropies for SU(A^) theory to up order 0(/3^) 
read 

sr-) . ^ (1±^) ^ 0(^y (22) 

sr-' ^ ^ (li^) + ^ Oin (23) 

^ (l±^ j ^ oin (24) 

As a result, the entanglement entropy of SU(A^) gauge theory in the strong coupling 
expansion is 


q{N=2) _ ^ 
‘^EE — ^2 

“± 

q{N=3) _ ^ 

‘^EE — ^2 

a_|_ 

q(N>3) _ ^ 
“^EE - 0 

LL I 


4iV^ N'‘ „,N‘ ^ 

-—+ 6^+0(^ 

4iV2 iV2 4 


+ 5S- 


{N=2) 




EE ’ 


A2 3 A3 

AN^ 

+ 3—+ 0(—) 


+ 5R 


iN=3) 


EE 1 


+ iS'F”- 


(25) 

(26) 
(27) 


A ' ^ A2 ^ A4 

where A± is the area of interface, a± is lattice spacing in the transverse space and where 
A = g‘^N is the t’Hooft coupling. As we argued below Eq. (12), this result is independent 
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of how the lattice is discretize, as long as the conical singularity is encircled by the same 
number of plaquettes in the action. The additional terms (5 S'ee are the positive contribution 
from the cosmological constant living on the d — 1 dimensional space transverse to the 
conical singularity in the Lagrangian, which we expect to cancel the negative term in the 


entanglement entropy, and will be explained in Sec. (IV). 

Likewise, the U{1) result can also be obtained using the same method; 

- O 


^EE 


4 3 ^, 1 , 

- 2 - 4 4 ” 

9 ^ 9 ^ 9 ^ 


+ 


( 28 ) 


Note that this is for the Yang-Mills 17(1) gauge helds, instead of the electromagnetic 17(1). 


IV. A 2-D COSMOLOGICAL CONSTANT COUNTERTERM 


In the result Eqs. (25 27) and (28), one hnds that the leading order 0(/3) term in the 


entanglement entropy has a negative contribution while all the subleading orders are positive. 
It turns out that, with the 2 dimensional conical structure in a 4 dimensional space time, 
we are allowed to introduce more local operators in the continuum action 


S= d‘^x±d‘^x\\ 


■-f^ + C4 + C2(5®(l||,T) 


( 29 ) 


The C 4 term is a four dimensional cosmological constant counter term which does not con¬ 
tributes to the entanglement entropy. However, the two dimensional cosmological constant 
counter term C 2 , living on the space transverse to the conical singularity and breaks the trans¬ 
lational symmetry on the cone, can contribute to the entanglement entropy |18( . Assuming 
C 2 is a smooth function of n, then 


C 2 = C 2 (n - 1) + O {{n- 1)^) , (30) 

where we C 2 = C 2 (/?) is a function of /3. We have made use of the fact that C 2 should vanish 
at n = 1 where translational symmetry is recovered. Therefore there is an extra unspecihed 
contribution to the entanglement entropy which also obey the area law: 

^47 = ( 31 ) 

We label the N dependence explicitly since different theories would have different C 2 counter 
terms. Also, is jS dependent. 
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The negative leading term in our entanglement entropy could in principle be compensated 


by the contribution from the two dimensional cosmological constant in (31). We remind the 
readers here that the negative leading term arises from the constant term in the lattice gauge 
field Lagrangian. Since the constant in the Lagragian for the central plaquettes is different 


from that for the non-central ones by a factor of 1/n, the effect is just like the C 2 in Eq.(29). 
If these constant terms are not included in the lattice Lagrangian, like what was done in 
some of the actions studied previously |T5l IT71120] , then the negative contribution to the 
entanglement entropy will not arise, and the Lagrangian can not be reduced to the usual 
Yang-Mills one at continuum limit. We argue that the C 2 and C 4 terms will always appear 
by renormalization even they are set to zero at certain renormalization scale. Setting them 
to be zero is equivalent to choosing specific values for these counter terms. Previously it 
was known that different choices of the boundary conditions gave different values for the 
entanglement entropy m- This corresponds to employing different regularization schemes, 
but those differences can be compensated by having different values for the counter terms 
for different regularization schemes used. 

We further speculate that the two-dimensional cosmological constant can show up in the 
entanglement entropy of scalar and spinor field theories as well. Also, they could play the 
role of a counter term to absorb the ultraviolet divergence of entanglement entropy and 
make entanglement entropy a finite physical quantity. 


V. COMPARISON WITH PREVIOUS RESULTS 


Our system is of infinite size in spatial and time dimensions. In na, the entanglement 
entropy of SU{N) gauge fields on the lattice of finite size at finite temperature is calculated in 
the Lagrangian formalism by means of group characteristic expansions. In 1+1 dimensions, 
when the periodic boundary condition is imposed on the spatial dimension, the entanglement 
entropy is given by 




ent. 




(32) 


where A is the total area of the 1+1 dimensional spacetime and a is the lattice spacing. At 


the limit of infinite area, the entanglement entropy reduces to 0, instead of Eqs. (25-27). 
While the free boundary condition is chosen in the spatial dimensions, the entanglement 
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entropy vanishes identically. We can obtain the same result if we set the conical singularity 
on a site so there is no central plaquette at all. So the leading order result for a system 
of periodic boundary condition in 1+1 dimensions comes from the conhguration with the 
plaquettes tiling the whole dimensions. But this conhguration gives vanishing entanglement 


entropy when the free spatial boundary condition is taken. Despite this, the result in Eq.(32) 
does not have the area law due to the nature of 1 spatial dimensional system, and hence 
does not have the expected form. 

Ref. [TU] considers the entanglement entropy of the ground states of the SU(2) Kogut- 
Susskind Hamiltonian |21| for the Wilson gauge theories at strong coupling limit, by in¬ 
cluding the edge states living on the boundary into the Hilbert space, such that the total 
entanglement entropy contains the contribution from the edge states. For d + 1 > 3 dimen¬ 
sions, the leading order entanglement entropy obtained by nm is 


‘S'ent. — — y (d — l)/3^(ln — + 1 + 21n2) 

CL LJ 


(33) 


where A± is the boundary area j2^. The entire leading contribution to the entropy in (33) is 
given by the entanglement entropy of the edge modes, which corresponds to the contribution 
of the consmological constant counterterm in our model. The non-local correlations of the 
d.o.f.’s in the two subregions are manifest only at higher order. Ref. jT2] also gives similar 
result: it is demonstrated via numerical simulation that, in the case of the Z 2 lattice gauge 
theory in three dimensions, the entanglement entropy is almost saturated by the entropy of 
the end points of the electric strings cut open by the boundary of the two subregions. In our 
language, their setup corresponds to locating the conical singularity on a site, so there is no 
central plaquette at all such that the leading contribution is coming from the counterterm. 


VI. SUMMARY 

To summarize, we have calculated the entanglement entropy of the SU{N) Yang-Mills 
gauge theories on the lattice under the strong coupling expansion in powers of /3 = 2N/g‘^. 
Using the replica method, our Lagrangian formalism maintains gauge invariance on the 
lattice. At 0(/3^) and 0{f3^), the entanglement entropy is solely contributed by the central 
plaquettes enclosing the conical singularity of the n-sheeted Riemann surface. The area law 
emerges naturally to the highest order 0(/3^) of our calculation. The leading 0(/3) term 
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is negative, which could in principle be canceled by taking into account the cosmological 
constant living in interface of the two entangled subregions. This unknown cosmological 
constant resembles the ambiguity of edge modes in the Hamiltonian formalism. We have 
further speculated that this unknown cosmological constant can show up in the entanglement 
entropy of scalar and spinor field theories as well. Furthermore, it could play the role 
of a counterterm to absorb the ultraviolet divergence of entanglement entropy and make 
entanglement entropy a finite physical quantity. 
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